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1  Introduction 

In  still  fluid,  heavy  particles  tend  to  fall  and  reach,  after  some  time,  a  terminal  velocity 
determined  by  the  balance  between  the  pull  of  gravity  and  the  viscosity  of  the  fluid.  Often, 
plankton  cells  are  slightly  heavier  than  sea  water,  and  as  a  consequence  they  tend  to  sink 
in  the  water  column.  This  sinking  could  then  take  them  below  the  euphotic  layer,  arresting 
photosynthesis. 

To  overcome  sinking,  plankton  have  devised  various  types  of  strategies.  One  is  swim¬ 
ming,  that  is,  an  active,  energy-consuming  self-propulsion  than  can  oppose  gravitational 
sinking.  Other  phytoplankters  can  modify  their  buoyancy,  becoming,  at  least  for  some 
time,  positively  buoyant.  Others  can  exploit  the  turbulence  in  the  mixed  layer,  using  it  to 
stay  suspended  for  longer  times.  This  latter  option  was  discussed  at  length  by  Margalef  [5], 
who  associated  the  level  of  turbulence  in  the  surface  layer  with  the  type  of  organisms  that 
are  environmentally  favoured.  The  terminal  velocity  of  spherical  impurities  with  the  same 
density  grows  with  the  squared  radius  of  the  falling  impurity.  Thus,  larger  plankton  cells 
tend  to  sink  faster  than  smaller  ones,  provided  they  can  be  approximated  by  a  somewhat 
spherical  shape.  When  the  turbulence  in  the  euphotic  layer  is  low,  Margalef  argued,  smaller, 
possibly  swimming  phytoplankters  such  as  flagellates  are  favoured  as  the  larger  species  sink 
more  rapidly  out  of  the  euphotic  zone.  On  the  other  hand,  when  the  level  of  turbulence 
is  high,  larger  plankton  can  benefit  from  turbulent  suspension  and  become  favoured  with 
respect  to  the  small  ones  (which  feel  diffusion  limitation  of  nutrients  more  severely). 

One  potential  problem  with  this  view  is  that  it  is  not  clear  whether  turbulence  really 
favours  a  prolonged  suspension  of  heavy  impurities.  For  example,  the  results  of  an  experi¬ 
ment  by  Ruiz  et  al  [8]  apparently  indicate  that  turbulence  make  heavy  particles  sink  faster 
than  in  still  fluid.  However,  other  experiments  considering  different  flow  configurations  gave 
the  opposite  result,  see  for  example  [1].  A  careful  study  of  this  issue  is  thus  required. 


2  Suspension  of  heavy  impurities 


The  non-dimensional  equation  of  motion  for  a  heavy  impurity  under  the  influence  of  gravity 
is  written  as 
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where  V  is  the  velocity  of  the  impurity,  u  is  the  fluid  velocity,  W  is  the  terminal  velocity  of 
the  impurity  in  still  fluid,  z  points  vertically  upwards  and  St  is  the  Stokes  number  defined 


as 


St 


2 

9h 


where  6  =  Pf/pp  is  the  ratio  between  the  density  of  the  fluid  and  the  density  of  the  impurity, 
L  is  a  typical  lengthscale  of  the  flow,  U  is  a  characteristic  flow  velocity  and  Re  the  Reynolds 
number.  The  terminal  velocity  W  is  found  by  equating  the  gravitational  force  to  the  Stokes 
drag,  and  its  non-dimensional  version  is  given  by 


2.1  Permanent  suspension  in  a  cellular  flow  field 

Stommel  [10]  studied  the  vertical  settling  of  small  heavy  impurities  in  a  simple  flow  configu¬ 
ration  composed  by  a  network  of  vortices  with  horizontal  axis;  all  the  dynamics  takes  plaxe 
in  the  vertical  plane  (x,  z)  and  the  vortices  are  infinitely  long  in  the  y  direction.  Stommel 
neglected  the  pressure  term  5Du/Dt  (which  we  shall  also  do  in  the  rest  of  this  section) 
and  neglected  also  the  inertia  of  the  impurities,  imposing  dV/dt  =  0  in  equation  (1),  which 
gives 

V  =  u  -  zW.  (3) 

Defining  the  two-dimensional  stationary  streamfunction 

T  =  2  (cos  x  +  cos  z) .  (4) 


The  flow  held  is  given  by 

u  =  (u,  w ) 

and  the  velocity  of  the  impurity  is 
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where  —  Wx.  The  solution  is  sketched  in  Figure  1.  The  how  is  characterized 

by  regions  with  closed  streamlines  of  \k,  surrounded  by  jet-like  channels  with  essentially 
downward  velocities,  owing  to  the  effects  of  gravity.  The  impurities  which  end  up  in  the 
downward  jets  fall  faster  than  in  still  fluid,  but  some  impurities  get  trapped  in  the  cells  with 
closed  streamlines  and  remain  suspended  forever.  Stommel  [10]  thus  argued  that  permanent 
suspension  is  possible  also  for  heavy  particles  and  that  turbulence  (or  cellular  how  helds) 
can  keep  heavy  impurities  afloat.  An  extension  of  this  approach  to  the  case  of  a  time- varying 
stream  function  was  discussed  in  [9]. 
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Figure  1:  The  motion  of  passively  advected  heavy  impurities  in  the  cellular  flow  field  con¬ 
sidered  by  Stommel  [10]. 


(a)  Circular  cell 


Figure  2:  (a)  Heavy  impurities  in  an  array  of  circular  flow  cells  drift  away  from  the  center 
owing  to  inertial  effects,  (b)  Impurities  in  a  (frozen)  flow  field  with  more  complex  streamlines 
can  stay  suspended  forever  (from  [7]). 


2.2  The  role  of  inertia 

According  to  Maxey  and  Corrsin  [6]  permanent  suspension  is  actually  not  possible  when 
particle  inertia  is  not  neglected.  Owing  to  the  centrifugal  force  acting  on  a  particle,  in  their 
motion  in  the  circulation  cells  the  impurities  will  spiral  out,  away  from  the  cell  center,  as 
illustrated  in  Figure  2a.  This  has  later  become  known  as  the  sling  effect:  heavy  particles 
tend  to  be  ejected  from  high  vorticity  regions  and  concentrate  in  high  strain  regions  be¬ 
tween  the  flow  cells;  closed  orbits  no  longer  exist  and  heavy  particles  eventually  end  up  in 
the  downdrafts  between  the  cells,  causing  an  increase  of  the  mean  settling  rate  compared  to 
what  happens  in  still  fluid,  in  keeping  with  the  results  of  the  experiment  mentioned  above  [8]. 

However,  consideration  of  flows  with  a  more  complex  streamline  pattern,  such  as  that 
depicted  in  figure  2b,  can  change  the  picture.  The  work  of  [7]  considered  a  spatially  ran¬ 
dom,  two-dimensional,  incompressible  stationary  flow  with  energy  spectrum  E(k)  oc  A;-'5/3. 
In  this  flow,  impurities  are  now  pushed  by  their  own  inertia  outwards  as  well  as  inwards, 
depending  on  the  streamline  curvature.  In  turn,  this  can  lead  to  the  emergence  of  closed 
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trajectories,  which  corresponds  to  a  limit  cycle  in  the  four-dimensional  impurity  phase  space 
(X,Y,Vx,Vy). 

Thus,  the  study  of  the  motion  of  small,  heavy,  spherical  particles  in  a  random,  station¬ 
ary  two-dimensional  flow  shows  that  permanent  particle  suspension  is  possible  when  the 
randomness  suitably  alters  the  curvature  of  the  streamlines  around  an  eddy.  A  stationary 
random  two-dimensional  flow  will  have  two  distinct  groups  of  impurities,  as  shown  in  figure 
3a.  The  particles  in  one  group  will  stay  very  close  to  where  they  started  (that  is,  they  stay 
suspended),  moving  on  limit  cycles  with  closed  trajectories.  The  impurities  in  the  other 
group  move  away  from  where  they  started,  reach  the  downdrafts  and  sink  with  an  average 
settling  speed  which  is  larger  than  their  terminal  velocity  in  still  fluid. 

Owing  to  inertial  effects,  identical  particles  can  thus  have  significantly  different  behav¬ 
ior  when  moving  in  a  random  steady  flow.  Some  particles  are  swept  into  the  downdrafts 
between  the  eddies  and  settle  at  rates  that  are  larger  than  in  still  fluid.  Other  particles 
remain  suspended,  moving  on  closed  trajectories,  so  that  their  vertical  motions  are  oscil¬ 
latory.  These  two  different  types  of  behavior  result  in  a  bimodal  distribution  of  vertical 
displacements.  As  a  result,  the  mean  settling  velocity  can  be  either  larger  or  smaller  than 
the  terminal  velocity  in  still  fluid  depending  on  the  fraction  of  suspended  particles. 

For  a  time- varying  random  two-dimensional  flow,  the  distinction  between  the  population 
of  suspended  impurities  and  that  of  rapidly  falling  ones  does  not  remain  so  strict  as  particles 
can  move  from  one  group  to  the  other.  In  this  case,  permanent  suspension  is  no  longer 
possible,  but  suspension  can  still  last  for  a  long  time.  Overall,  these  results  confirm  the  view 
of  Stommel  and  Margalef,  who  claimed  that  turbulence  could  indeed  help  heavy  impurities 
to  stay  suspended  for  longer  times  than  in  still  fluid. 

3  Phytoplankton  competition  in  turbulent  waters 

As  turbulence  affects  the  sinking  rate  of  phytoplankton,  one  may  go  back  to  the  old  Mar¬ 
galef  question  and  ask  how  populations  of  phytoplankton  species  with  different  size  respond 
to  a  turbulent  environment.  Larger  plankton  species  seem  to  be  favoured  in  conditions  of 
strong  turbulent  suspension  while  smaller  phytoplankton  can  thrive  in  quieter  conditions. 

To  explore  this  issue,  we  follow  [3]  and  consider  two  phytoplankton  species  with  radius 
ai  and  a2  that  compete  for  the  same  nutrient,  N.  The  population  dynamics  is  described 
by  the  following  (non-dimensional)  NP  system: 

DN  N  Pt  NPo 

=1  -  — — Q(ai)  -  f3(a2)  —  ~NQ(a2)  +  n ntu  [PiQ{a{)  +  P2Q{a2)\ 

=^ai^k  +  N  ~  +  s(ai^  Pl  (5) 

DP2  qi  .  np2 

~W  =l3MkTN  - [m  +  5(02)1  P2' 
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AY(Lo) 

(a)  Frozen  flow  field 


AY(Lo) 

(b)  Time-evolving  flow  field 


Figure  3:  (a)  In  a  frozen  turbulent  flow  some  impurities  can  stay  suspended  forever,  while 
others  move  away  from  their  initial  position  and  fall  faster  than  in  still  fluid.  As  a  result, 
the  distribution  of  particle  displacements  is  bimodal.  The  three  curves  show  the  probability 
distributions  of  the  vertical  displacement  of  falling  impurities  at  different  times  and  the 
arrows  indicate  the  displacement  that  the  impurities  would  have  had  if  they  moved  at  the 
Stokes  settling  velocity  W.  On  the  left  of  the  plot,  there  are  particles  which  stay  suspended 
forever  and  their  displacement  from  the  initial  position  does  not  grow  with  time,  (b)  In 
a  time-evolving  random  field  there  is  no  strict  separation  between  suspended  and  falling 
particles  and  permanent  suspension  is  not  possible,  but  a  large  fraction  of  the  impurities 
can  fall  more  slowly  than  in  still  fluid.  From  [7]. 


where  Pi  and  P2  are  the  concentrations  of  the  two  species  of  phytoplankton,  f3  is  a  growth 
rate,  k  is  a  limiting  term  for  phytoplankton  growth,  m  is  a  size-independent  natural  phy¬ 
toplankton  mortality  rate,  s  is  the  size-dependent  sinking  rate,  and  Q  is  the  nutrient  quota 
of  phytoplankton,  used  to  transform  biomass  into  nutrient  amount. 

The  terms  in  the  nutrient  equation  of  system  (5)  describe  two  sources  of  nutrients:  a 
constant  external  nutrient  source,  /,  and  the  nutrients  regenerated  from  dead  phytoplank¬ 
ton,  as  described  by  the  recycling  term  [Pi  Q(a\  )  +  FbQ^)]  with  //jv  <  1-  Nutrients 
are  consumed  by  phytoplankton  at  a  rate  (3 ,  describing  the  phytoplankton  efficiency  to  up¬ 
take  nutrient  which  is  assumed  to  be  a  function  of  phytoplankton  size. 


In  their  study,  Jiang  et  al  [3]  recalled  that  there  is  a  complex  relationship  between  the 
efficiency/growth  rate  /?  and  the  size  of  phytoplankton  species.  Observations  show  that  the 
growth  rate  of  relatively  large  phytoplankton  species  decreases  with  increasing  size  while 
the  opposite  tends  to  be  true  for  small  phytoplankton  species.  This  suggests  that  there 
should  be  an  optimal  phytoplankton  size  at  which  the  growth  rate  is  maximized.  A  simple 
analytic  form  of  the  phytoplankton  growth  rate  that  has  these  characteristics  is 


/3(a) 


a 

a\ a2  +  0:2a  +  0:3 


(6) 


where  aq,  0:2,  and  <33  are  positive  constants.  The  form  of  the  growth  rate  described  in  (6) 
is  only  meant  to  capture  the  qualitative  structure  of  observed  phytoplankton  growth  rates 
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Figure  4:  Outcome  of  species  competition  in  a  turbulent  environment,  parameterized  by 
the  values  of  the  exponent  a  and  the  coefficient  so-  Red  (in  the  upper  left  area)  indicates 
dominance  of  the  larger  species  with  a  =  1.4,  blue  indicates  dominance  of  the  smaller  species 
with  a  =  0.7.  Cell  radii  are  in  non-dimensional  units. 

and  is  not  necessarily  universal. 

The  growth  of  phytoplankton  populations,  thanks  to  nutrient  consumption,  is  balanced 
by  plankton  mortality,  m,  and  plankton  sinking,  s(a).  In  the  phytoplankton  competition 
model  (5),  the  mortality  rate  rn  is  assumed  to  be  independent  of  phytoplankton  size,  un¬ 
like  the  sinking  term.  In  the  absence  of  turbulence  and  assuming  spherical  particles,  the 
settling  velocity  of  an  individual  phytoplankter  can  be  derived  from  Stokes  law  obtaining 
s(a )  =  so®2;  this  form  was  assumed  in  [3]  with  so  =  2.  With  this  choice,  it  was  found  that 
the  most  favoured  species  is  always  the  smaller  one,  and  that  there  should  be  a  natural 
tendency  towards  dominance  of  smaller  phytoplankton  species.  In  this  view,  the  higher  ef¬ 
ficiency  in  nutrient  uptake  of  the  larger  species  is  not  enough  to  balance  their  faster  sinking 
rate,  leading  to  a  dominance  of  the  smaller  organisms. 

In  reality,  the  sinking  rate  depends  on  the  level  of  turbulence,  which  can  alter  the  above 
result.  In  a  subsequent  study,  Demonte  et  al  [2]  kept  all  parameter  values  as  in  [3]  but 
assumed  the  sinking  rate  to  be  proportional  to  a  power  of  the  phytoplankton  cell  radius,  o, 
through  a  simple  relationship 


s(a)  =  sqo'7 


(7) 


where  so  is  the  coefficient  and  a  is  the  exponent.  The  motivation  of  this  choice  is  that 
turbulence  can  alter  the  simple  dependence  of  sinking  velocity  on  radius  observed  in  still 
fluid.  What  is  found  in  this  case  is  that  the  outcome  of  species  competition  depends 
on  the  properties  of  turbulence,  that  is,  on  the  assumed  values  of  sq  and  a.  Figure  4 
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shows  the  outcome  of  competition  for  two  species  with  (non-dimensional)  radius  0.7  and 
1.4,  respectively.  For  the  parameter  values  chosen  in  [3],  the  smaller  species  is  favoured. 
However,  by  varying  the  values  of  a  and/or  so,  one  can  get  the  opposite  result.  Even  by 
keeping  a  =  2  but  varying  the  coefficient  so,  one  can  obtain  dominance  of  the  larger  species. 
This  is  consistent  with  the  view  that  the  outcome  of  competition  of  phytoplankton  species 
with  different  size  depends  on  the  properties  of  the  turbulent  environment,  and  it  seems  to 
further  confirm  Margalef’s  view. 

4  Vertical  velocities  in  mesoscale  flows 

In  the  previous  sections  we  have  seen  that  turbulent  suspension  can  be  an  important  mech¬ 
anism  in  plankton  dynamics.  We  can  ask  whether  these  effects  are  visible  also  on  scales 
larger  than  those  of  homogeneous  and  isotropic  turbulence:  for  example,  whether  they  can 
play  a  role  at  the  ocean  mesoscales  where  fronts,  vortices  and  waves  dominate  the  flow,  as 
discussed  in  the  first  lecture. 

Although  a  final  answer  is  still  missing,  an  important  point  concerns  the  structure  and 
intensity  of  the  vertical  velocity  field  at  the  mesoscales.  If  there  are  strong  and  organized 
vertical  velocities,  then  these  can  also  induce  suspension  (in  upwelling  regions)  and  rapid 
sinking  (in  downwelling  areas),  thus  affecting  plankton  dynamics,  besides  their  role  in  de¬ 
termining  the  nutrient  input  to  the  euphotic  layer. 

In  proximity  of  fronts,  vertical  velocities  can  be  very  intense,  creating  a  complex  pattern 
of  upwelling  and  downwelling.  However,  other  mesoscale  structures  can  also  play  a  role:  co¬ 
herent  vortices,  in  particular,  have  been  shown  to  be  associated  with  a  complicated  pattern 
of  strong  vertical  velocities  [4].  In  their  study,  Koszalka  et  al.  have  considered  a  primitive 
equation  model,  using  the  ROMS  ocean  model,  with  doubly  periodic  boundary  conditions 
and  idealized  surface  forcing.  The  simulations  show  the  emergence  of  coherent  vortices, 
similarly  to  what  happens  in  two-dimensional  or  quasi-geostrophic  (QG)  turbulence.  The 
horizontal  dynamics  and  transport  properties  of  this  flow  are  not  too  dissimilar  from  what 
happens  in  QG  turbulence,  suggesting  that  approaches  based  on  simple  QG  models  capture 
a  relevant  portion  of  the  horizontal  dynamics.  However,  the  vertical  velocities  associated 
with  the  vortices  in  the  primitive  equation  setting  can  be  very  strong,  at  variance  with  the 
situation  found  for  QG  vortices. 

Figure  5a  shows  a  slice  of  the  vorticity  field  in  one  of  the  primitive  equation  simulations 
discussed  in  [4],  and  figure  5b  shows  the  vertical  velocity  field  associated  with  one  of  the 
vortices.  The  vertical  valocities  in  proximity  of  the  vortex  are  very  intense,  and  can  lead 
to  various  consequences  for  the  marine  ecosystem,  ranging  from  enahnced  nutrient  input 
to  the  suspension  of  phytoplankters  in  the  upwelling  areas  and  their  rapid  sinking  in  the 
downwelling  regions. 
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Figure  5:  Panel  (a):  The  surface  vorticity  field  produced  by  a  numerical  simulation  of  forced- 
dissipated  primitive  equations  with  doubly  periodic  boundary  conditions  in  a  domain  with 
lateral  size  256  km  and  depth  1000  meters.  Forcing  is  provided  by  idealized  surface  (wind) 
forcing  and  relaxation  to  a  given  density  profile.  The  relative  vorticity,  C,  is  normalized  by 
the  Coriolis  parameter,  /.  Panel  (b):  The  vertical  velocity  field,  in  m/day,  at  a  depth  of  78 
meters,  for  the  vortex  indicated  by  the  white  segment  in  panel  (a).  From  [4], 

5  Conclusion 

In  these  lectures  I  tried  to  illustrate,  by  various  examples,  how  turbulence  affects  aquatic 
organisms.  The  overall  picture  shows  that  turbulence  is  an  important  component  of  these 
ecosystems  and  it  must  be  properly  considered  to  obtain  a  quantitative  description  of  ecosys¬ 
tem  functioning.  From  the  small  scales  of  homogeneous  and  isotropic  turbulence,  where 
individual  plankton  cells  swim,  sink  and  are  transported  by  the  fluid  flow,  to  the  ocean 
mesoscales  where  vortices,  fronts  and  waves,  with  their  associated  upwelling  and  down- 
welling  regions  determine  the  plankton  distribution,  to  the  large  scales  where  the  basin 
circulation  controls  the  overall  structure  of  the  marine  ecosystem,  fluid  dynamics  is  at  the 
core  of  the  biological  processes  occurring  in  aquatic  environments.  There  are  still  many 
open  questions,  some  of  which  can  possibly  be  addressed  by  the  methods  of  GFD  and  of 
dynamical  system  analysis  as  illustrated  here. 
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